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Calculations  of  surface  wave  and  radar  cross  section  modulation  induced  by  a  spatially  varying  surface 
current  have  been  compared  with  field  measurements  made  with  SAR  ( X -  and  L-  bands),  real  aperture 
radar  (Y-band),  laser  slope  gauge,  and  CCD  video  camera  during  the  Joint  Canada/U.S.  Ocean  Wave 
Investigation  Project  (JOWIP)  which  took  place  in  August  1983.  The  comparison  reveals  that  (i) 
many  existing  wind-relaxation  models  underpredict  the  hydrodynamic  effect  of  the  current,  (ii)  simple 
Bragg  scattering  model  underpredicts  radar  backscatter  modulation  for  higher  frequency  radars  (e.g., 
A"- band)  but  appears  acceptable  at  lower  radar  frequencies  (e.g.,  I- band),  (iii)  the  discrepancy  between 
measurements  and  calculations  at  high  Tadar  frequency  is  reduced,  but  not  eliminated,  when  the  effects 
of  long  surface  waves  have  been  accounted  for. 
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1.  INTRODUCTION 

One  of  the  major  objectives  in  the  Joint  Canada/U.S.  Ocean  Wave  Investigation  Project 
(JOWIP)  was  the  evaluation  of  the  capabilitites  and  limitations  of  theoretical  models  for  predicting  the 
modulation  of  real  and  synthetic  aperture  radar  backscatter  from  the  ocean  by  a  variable  surface 
current  induced  by  internal  waves,  flow  over  bottom  topography,  or  other  wavelike  or  non-wavelike 
disturbance. 

JOWIP,  or  the  Georgia  Strait  experiment  as  it  is  commonly  known  in  the  community,  has  the 
potential  of  making  a  substantial  contribution  to  the  understanding  of  ocean  physics  and  its  remote 
sensing  by  synthetic  aperture  radar  (SAR)  in  that  the  X-  and  £-band  simultaneous  airborne  SAR 
measurements  are  supported  by  a  full  set  of  ground  truth  measurements,  providing  in  situ 
measurements  of  subsurface  characteristics  (CDT  profiles),  surface  characteristics  (x-y  current,  wind 
velocity,  surface  wave  slope),  real  aperture  radar  returns  (X-band  radar  supported  by  optical  specular 
detector),  optical  returns  (CCD  video  camera),  stereo  photographs,  and  airborne  photographs.  The 
SAR  overflights  and  ship  crossings  are  co-ordinated,  synchronized  and  tracked. 

The  wealth  of  quantitative  ground  truth  data  taken  during  SAR  overflights  makes  JOWIP 
uniquely  qualified  to  validate  and  assess  a  theoretical  model.  Measurements  of  surface  current,  surface 
wave  slope,  and  real  aperture  radar  (RAR)  backscatter  provide  all  of  the  data  needed  to  test  each  of 
the  major  components  in  such  a  model.  The  JOWIP  data  collection,  unlike  for  example,  SARSEX 
data  which  lacks  direct  surface  current  measurements  and  real  aperture  radar  data,  can  be  used  to  test 
the  hydrodynamic  model  by  comparing  with  the  current  and  wave  slope  data,  the  radar  scattering 
model  by  comparing  with  the  wave  slope  and  RAR  backscatter  data,  and  the  radar  imaging  model  by 
comparing  with  the  wave  slope,  RAR  and  SAR  data.  Without  these  intermediate  tests  of  the 
components,  the  large  number  of  variables  in  an  end-to-end  SAR  simulation  model  makes  it  all  too 
easy  to  improperly  attribute  successes  or  failures  of  one  component  to  another,  or  to  arrive  at 
misleading  conclusions  due  to  cancelling  effects  of  simplifications  made  within  different  components. 

In  this  paper,  we  test  the  models  for  (a)  hydrodynamic  interaction  of  a  surface  wind  wave  field 
with  surface  currents,  (b)  radar  scattering  from  the  modulated  surface  waves.  We  compare  results  of 
model  calculations  obtained  with  the  TRW  Hydrodynamic  Interaction  Code  (TRW-HIC)  to  remote 
measurements  obtained  by  synthetic  aperture  radar  (SAR)  and  to  surface  measurements  obtained  by 
ship-borne  slope  gauge  and  real  aperture  CW  radar  (RAR).  The  comparison  focuses  on  the  validity  of 
several  important  elements  of  the  theoretical  model: 
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a.  The  model  for  wind  effects  (relaxation  model); 

b.  The  model  for  the  equilibrium  spectrum  in  the  absence  of  a  variable  surface 
current; 

c.  The  effects  of  certain  simplifying  assumptions  in  the  model; 

d.  The  Bragg  model  for  radar  baclucatter; 

e.  The  long  wave  field  effects. 


Computations  were  made  for  a  number  of  runs  to  address  the  above  issues.  Two  cases 
(JOWIP  1-1  and  TRW  8 IB),  were  selected  for  detailed  quantitative  comparisons. 

Shortly  after  the  completion  of  the  work  reported  here,  a  new  radar  ocean  imaging  model  was 
proposed  by  Holliday  ti  al  (1986).  Comparison  between  the  results  presented  here  and  results  obtained 
by  using  this  new  model  will  be  the  subject  of  a  following  paper. 


2.  MODEL  FOR  BRAGG  WAVE  MODULATION  BY  A  VARIABLE  CURRENT 

The  Bragg  model  proposes  that  the  baclucatter  of  electromagnetic  radiation  from  the  ocean 
surface  is  mainly  due  to  waves  of  a  specific  wavelength  through  a  grating  like  mechanism.  For  a  radar 
with  wave  number  kT  and  incidence  angle  6,  the  model  predicts  that  only  water  waves  of  wavenumber 
kg=2kTsinO  in  the  line  of  sight  contribute  to  the  baclucatter.  The  radar  cross  section  <r  is  predicted  by 
the  Bragg  model  to  be  proportional  to  F(kg),  where  F  is  the  spectral  density  function  of  water  waves 
(Wright,  1966).  The  model  assumes  a  spatially  homogeneous  wave  field,  and  for  sufficiently  steady 
conditions,  F  is  given  by  the  equilibrium  wave-height  distribution  Fef(k)  for  the  prevailing  wind  (see, 
for  example,  Pierson  and  Stacy,  1973). 


A  straightforward  extension  of  the  model  to  the  case  where  the  surface  wave  spectrum  is 
spatially  modulated  by  a  variable  surface  current  U(x,t)  allows  the  identification  of  relative  spatial 
changes  in  baclucatter  cross  section  with  relative  spatial  changes  in  the  spectral  content  at  the  Bragg 
resonant  wavenumber.  If  a  dimensionless  spectral  density  perturbation  function  f(k,x,i)  associated 
with  U(x,t)  is  defined  as 


(1) 


/(*,m 


F(k,l,t)  -  Feq(j) 
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presence  of  the  current  is 

(2)  [i  +  /(*«.*.<)]. 

where  a  0(kB)  is  the  steady  state  Bragg  backscatter  cross  section  in  the  absence  of  a  surface  current. 
The  relative  modulation  of  the  backscatter  cross  section  is  then  the  nondimension&l  perturbation  in  the 
surface  wave  spectrum  of  the  Bragg  wave  component: 

f  =  lCk‘x,)- 

**  o(*B' 


ts 

(3) 

m  4 

where  Sw 

V 

(4) 

£ 

V' 

(5) 

where 

A.  Governing  Equations. 

In  the  absence  of  source  terms  (including  wind  energy  input  and  dissipation),  the  local  wave 
action  N(k,z,t)  is  conserved.  In  the  presence  of  wind,  N  obeys  the  equation: 

dN  _  ON  ,  A  v  iz  ,  a  „  dk  c 
—  =  -S7  +  +  vkN-t Ti  =  sw, 


it  di 


dz  _  dw 

di  ~  dk 


dt 


dk 

dt 


— 
di  ’ 


(6) 


u>  =  u>'  +  k-U(z,t) 


is  the  local  dispersion  relation,  and  u1  is  the  intrinsic  frequency  given  by 

(7)  u/J  =  gk(l+Tk3), 

with  T  the  surface  tension  coefficient,  and  k  =  |i|. 

The  action  density  N(k,z,t)  is  related  to  the  spectral  density  function 
F(k,z,i)  by  the  relation 
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(8) 


N(U,t)  =  j  F(k,x,t). 


In  the  absence  of  a  current  pattern,  the  surface  wave  field  is  assumed  to  be  established  by  a 
constant  wind  blowing  for  a  long  time,  and  is  described  by  a  spatially  homogeneous  equilibrium 
spectral  density  function  fef(k),  which  is  determined  heuristically. 

The  wind  source  term  Sw  is  assumed  to  take  the  form 


0) 


sw  =  p  j  Fli~  iw,r]. 


where  a  can  be  1  or  2  depending  on  whether  the  restoring  force  is  quadratic  or  cubic.  This  source  term 
acts  to  (a)  develop  a  flat  surface  into  the  equilibrium  spectrum  under  a  given  wind,  and  (b)  restore  any 
departure  from  this  equilibrium,  at  the  same  rate;  this  “relaxation”  rate  0  being  a  function  of  wind 
and  wavelength.  The  relaxation  rate  0  is  also  determined  heuristically,  based  mainly  on  wave  growth 
data  in  the  laboratory  and  the  ocean  (cf.  Section  2Bb). 

The  surface  current  U(x,t)  is  taken  to  be  a  steady  pattern  propagating  in  the  x  direction  with 
a  speed  V.  In  the  frame  of  reference  moving  with  V,  the  current  pattern  is  described  by  U(i‘)  = 
U(z—Vt,y),  and  the  resultant  spatially  inhomogeneous  equilibrium  spectral  density  function  can  be 
written  in  terms  of  Fef  and  of  the  dimensionless  spectral  perturbation  /  defined  in  Equation  (1)  and 
evaluated  at  (k, I*).  (Henceforth,  we  will  drop  the  *  when  the  context  is  clear). 

Incorporating  Equations  (1),  (8)  and  (9)  into  Equation  (3)  yields  the  evolution  equation  for 

/(M) 


(10) 


Tt  = 


w(  “+2r')  + 


1  dk  dkfkdu 
k  di  'di  \»  dk 


F «♦  9k  ) 


According  to  the  Bragg  scattering  model,  the  solution  of  Equation  (10)  along  the  wave 
trajectories  defined  by  Equations  (4)  and  (5),  gives  the  relative  changes  in  backscatter  cross  section  for 
a  given  wind  and  a  given  surface  current  distribution. 


B.  Modeling  issues. 


The  approach  described  above  requires  that  the  current  pattern  propagation  and  the  wind 
conditions  be  reasonably  steady.  Assuming  this  to  be  fulfilled  in  the  cases  where  the  model  is  applied, 
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there  are  several  issues  that  are  relevant  to  the  validity  of  the  model. 


Model  for  Equilibrium  Spectrum. 


In  the  case  of  the  Georgia  Strait  experiments,  preliminary  comparison  with  background  slope 
gauge  measurements  indicated  that  the  Pierson-Stacy  (1973)  model  is  in  reasonable  agreement  with 
data  in  the  energy-containing  part  of  the  spectrum.  Based  on  that  comparison,  this  equilibrium 
spectrum  was  used  in  the  calculations  presented  here.  Detailed  comparisons  in  the  high  wavenumber 
range  and  on  directional  dependence  have  not  been  made. 


b.  Model  for  the  Wind  Source  Term  and  Relaxation  Rate. 


There  are  several  primarily  heuristic  models  available  for  the  relaxation  rate  0.  They  are 
based  mainly  on  different  sets  of  wave  growth  data  in  the  laboratory  and  in  the  ocean,  with  some 
overlapping  in  the  parameter  regime.  Models  that  we  have  implemented  are: 


a.  Inoue  (1966)  model, 

b.  Miles  (1959)  model, 

c.  Snyder  and  Cox  (1966)  model, 

d.  Snyder  et  al  (1981)  model, 

e.  Mitsuyasu  and  Honda  (1982)  model, 

f.  Hughes  (1978)  model, 

g.  Plant  (1982)  model. 


The  first  five  models  listed  are  to  be  used  with  a  cubic  restoring  force  in  the  wind  source  term 
(i.e.,  a=2  in  Equation  (9)),  while  the  last  two  models  are  for  a  quadratic  (a=l)  restoring  force.  In  all 
of  these  models,  for  a  given  wind,  the  shorter  waves  have  a  larger  relaxation  rate  (shorter  relaxation 
time).  Hence,  they  are  more  responsive  to  wind  action.  However,  the  actual  values  for  0  are  spread 
over  an  order  of  magnitude  among  the  different  models,  as  illustrated  in  Figure  l.a  for  the  case  of  a 
5  m/sec  wind.  Also,  there  is  presently  no  reliable  data  to  support  any  specific  directional  dependence 


For  a  given  wind  speed,  the  Inoue  model  gives  the  largest  relaxation  rate,  thus  the  strongest 
wind  effect.  When  compared  with  the  Snyder  ei  al  (1981)  model  (labelled  “Sea”  in  Figure  l.a),  Inoue’s 
model  has  a  relaxation  rate  about  one  order  of  magnitude  larger,  and  hence  predicts  a  modulation 
about  one  order  of  magnitude  smaller.  A  sample  of  results  obtained  using  various  0  models  is  shown 
in  Figure  l.b. 
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Preliminary  comparison  with  Georgia  Strait  data  clearly  indicated  that  calculations  using  the 
Inoue  model  were  underpredicting  the  measurements.  For  this  reason,  the  model  yielding  the  largest 
perturbation  -  which  is  the  Snyder  et  al  (1981)  model  -  was  selected  over  the  Inoue  (1967)  model. 
This  selection  has  been  confirmed  by  detailed  comparison.  The  Hughes  (1978)  model,  when  used  in 
conjunction  with  a  quadratic  restoring  force  term  (which  effectively  reduces  the  relaxation  rate  by  a 
factor  of  2  for  small  departures  from  equilibrium),  yields  results  comparable  to  the  Snyder  et  al  model 
for  small  modulations. 


When  applied  to  the  prevailing  conditions  of  these  experiments,  the  various  0  models  would 
have  produced  results  spread  over  a  range  larger  than  one  order  of  magnitude.  The  selection  of  the 
Snyder  et  al  (and  Hughes)  models  over  the  Inoue  model  is  a  very  important  conclusion  from  the 
JOWIP  data.  It  implies  that  previous  calculations  which  used  the  Inoue  model  may  have 
underpredicted  surface  modulations  induced  by  internal  waves  by  up  to  an  order  of  magnitude. 
However,  it  should  be  pointed  out  that  this,  as  many  of  the  other  conclusions,  is  based  on  a  limited 
number  of  cases  analysed.  More  data  (including  JOWIP  data  not  yet  processed)  should  be  analysed  to 
fortify  or  challenge  a  conclusion  as  important  as  this. 


c.  Simplifying  Assumptions. 

Several  popular  approaches  to  the  solution  of  Equation  (10)  involve  further  simplifications. 
We  make  a  few  observations  about  the  effect  of  some  of  these  simplifications  in  the  application  to  the 
Georgia  Strait  experiments. 

i)  Linearization  about  the  Equilibrium  Spectrum. 

When  the  current  and  its  variations  are  small,  the  induced  perturbations  are  also  expected  to 
be  small.  The  a  priori  assumption  that  f  «  1  justifies  linearization  of  the  equation  about  the 
equilibrium  spectrum.  Instead  of  solving  Equation  (10)  which  is  nonlinear  in  /,  the  governing  equation 
becomes: 


(ID 


d£ 

di 


<*0f  + 


du  k  du  .  k  dF"~ 
dx  [_  w  dk  Ft1  dk 


Comparison  between  the  full  and  linearized  calculations  for  surface  currents  corresponding  to 
the  Georgia  Strait  experiment  show  good  agreement  for  those  portions  of  the  record  where  current 
gradients  ore  very  weak.  However,  where  the  modulation  /  is  not  very  small,  say  larger  than  10%,  the 
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i _ . _ i _ i 

solution  to  the  linearized  equation  departs  considerably  from  that  of  the  full  equation,  leading  in  many 
cases  to  unphysical  negative  values  for  the  spectral  density  F  (i.e.,  /  <  -1). 

ii)  Linearised  Characteristics. 

Another  simplification  is  to  employ  straight  characteristics  rather  than  curved  characteristics 
when  integrating  the  spectral  equation.  The  validity  of  this  assumption  can  be  pre-tested  by  plotting 
the  characteristics  before  performing  the  calculations.  An  example  of  such  a  plot  for  JOWIP  1-1  is 
shown  in  Figure  2.  It  can  be  seen  that  in  certain  wavenumber  ranges,  the  waves  are  actually  “trapped” 
by  the  current  pattern.  For  both  L-  and  A- band,  the  characteristics  are  strongly  curved,  indicating 
that  the  straight-characteristics  assumption  fails  for  this  case. 

iii)  Beta-Dominant  Approximation. 

The  so-called  “beta-dominant”  approximation  assumes  that  the  wind  restoration  mechanism  is 
dominant  in  the  sense  that  it  restores  the  wave  spectrum  to  a  new  equilibrium  instantaneously.  This 
assumption  leads  (Alpers  and  Hennings,  1984)  to 


(12) 


15  M 

a0  dz 


The  simplicity  of  the  expression  makes  this  a  popular  means  for  obtaining  a  “first  estimate”  on 
the  modulations!  magnitude. 


Since  the  wind  relaxation  rate  0  increases  with  decreasing  wavelength,  this  approximation  is 
expected  to  become  applicable  as  the  waves  of  interest  become  appropriately  short.  This  is  supported 
by  comparison  with  results  obtained  with  the  full  equations.  For  example,  in  the  case  of  JOWIP  1-1, 
the  beta-dominant  approximation  failB  at  L-band,  yielding  results  3  times  larger  than  those  from  the 
full  equations,  and  producing  negative  spectral  density  values.  For  A-band,  the  beta  dominant 
approximation  is  acceptable,  yielding  results  about  20%  larger  than  those  obtained  from  the  full 
equations. 


If  the  Inoue  relaxation  model  (which  has  a  0  value  one  order  of  magnitude  larger)  were  used, 
then  the  beta-dominant  solutions  approximate  the  solutions  of  full  equations  well  for  this  data  set  — 
both  being  about  an  order  of  magnitude  smaller  than  the  results  using  the  Snyder  et  al  model.  This 
suggests  that  the  apparent  success  attributed  to  the  beta-dominant  approximation  in  the  past  was  the 
result  of  adopting  a  large  /3-value  model.  With  the  selection  of  the  Snyder  et  al  model  over  the  Inoue 
model,  the  acceptability  of  the  beta-dominant  approximation  has  been  greatly  reduced. 
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C.  TRW-HIC  Implementation. 

Our  implementation  allows  the  user  to  specify  the  desired  equilibrium  wind  wave  spectrum,  to 
choose  among  any  of  the  relaxation  rate  models  listed  in  Section  2Ab,  to  use  either  a  quadratic  or  a 
cubic  wind  restoration  term,  and  to  impose  any  of  the  additional  simplifying  assumptions  listed  in 
Section  2Ac. 

The  main  output  is  the  dimensionless  spectral  perturbation  from  which  the  full 

spectrum  can  be  reconstructed  and  root-mean-square  wave  height  and  wave  slope  obtained. 
Additionally,  these  hydrodynamic  disturbances  can  then  be  converted  into  sensor  responses:  radar 
cross-section  modulation,  surface  glitter  pattern,  surface  imagery  of  sky  reflection,  and  geometric  effects 
on  radiometric  returns. 

For  the  results  presented  in  this  paper,  the  Pierson-Stacy  spectrum  was  used.  Also,  we  have 
chosen  the  wind  restoration  term  to  be  cubic  [«.e.,  a=2  in  Equation  (9)].  Given  the  evidence  referred 
to  in  Section  2Bb,  the  relaxation  rate  was  chosen  to  be  ^s  given  by  Snyder  ti  al  (1981).  For 
consistency,  0's  angular  dependence  was  taken  to  be  the  same  as  in  the  Pierson-Stacy  spectrum. 

None  of  the  simplifying  assumptions  of  Section  2Bc  was  exercised  in  the  generation  of  the 
results  presented  here.  However,  due  to  the  favorable  conditions  of  the  cases  chosen,  the  unidirectional 
version  of  the  code  was  used. 


3.  APPLICATION  TO  GEORGIA  STRAIT  EXPERIMENTS. 


A.  Geometric  Considerations. 

The  DREP  current  meters,  the  DREP  slope  gauge,  the  APL  video  camera  and  the  TRW  X- 
band  CW  RAR  are  all  ship  mounted.  These  measurements  comprise  the  set  of  ground-truth  data 
along  a  trajectory  traversed  by  the  ship.  These  ship- borne  instruments,  however,  do  not  measure  the 
exact  same  spot  of  water  surface. 

The  ERIM  SAR  provides  an  area  image.  A  cut  in  the  SAR  image,  colocated  with  the  ship 
trajectory  to  the  extent  possible,  is  used  for  comparisons  with  other  measurements  and  with 
calculations. 

The  calculated  results  of  spectral  perturbation  at  the  Bragg  resonant  wavenumber  are 
compared  with  the  slope  gauge  output  bandpassed  around  the  Bragg  component,  the  X-band  real 
aperture  result  (when  available),  the  video  data  (when  available)  and  the  SAR  cut.  To  the  extent  that 
these  various  output  measurements  are  not  spatially  or  temporally  colocated,  there  is  some  intrinsic 
uncertainty  in  the  comparison  itself  which  must  be  borne  in  mind  by  the  reader. 

The  input  quantities  for  a  calculation  are: 

a.  Wind  speed  and  direction  (from  meteorological  package), 

b.  Current  magnitude  and  direction  (from  current  meter), 

c.  Current  pattern  propagation  speed  and  direction  (both  from  delays  between  bow 
and  stern  current  meters  and  from  consecutive  SAR  images). 


These  vectors  are  resolved  in  the  direction  of  the  ship  motion.  The  calculations  are  therefore 
unidirectional  along  the  ship  track. 
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B.  Quantitative  Comparison  with  Data. 

Baaed  on  a  quick-look  analysis,  several  runs  were  selected  for  detailed  data  processing, 
including  JOWIP  1-1,  12-3/4,  14-4,  15-7  and  TRW  81B,  all  of  which  are  for  currents  induced  by 
naturally  occurring  internal  waves.  Among  the  most  important  selection  criteria  were  the  availability 
of  good  data  and  the  distinctive  characteristic  of  the  internal  wave  signals.  Calculations  were 
performed  for  all  these  cases.  Quantitative  comparisons  with  measurements  were  made  to  the  extent 
possible.  In  the  following,  detailed  comparisons  between  calculations  and  data  are  shown  for  JOWIP 
1-1  and  TRW  81B  only;  however,  the  general  conclusions  are  supported  by  the  other  cases  examined. 


a.  JOWIP  RUN  1-1. 

This  run  corresponds  to  data  taken  between  17:44:31  and  18:26:19  PDT  on  July  22,  1983.  In 
this  case,  the  wind  speed  was  recorded  at  10  m  height  as  3.4  m/s  blowing  toward  95'T.  The  ship  was 
heading  98*T  at  an  approximately  constant  speed  of  1  m/s.  The  ship  crossed  a  pattern  of  internal 
waves  associated  with  a  shallow  and  sharp  thermocline  with  a  maximum  Brunt-Vaisala  frequency  of 
60  cycle  per  hour  at  a  depth  of  2.5  m.  The  surface  current  as  a  function  of  distance  along  the  track  is 
shown  in  Figure  2a.  The  propagation  direction  of  the  pattern  was  determined  from  the  SAR  images  to 
be  about  1 10’T,  and  the  speed  from  the  time  delay  between  bow  and  stem  current  meters  to  be  0.01 
m/s.  (Velocities  are  all  land  referenced.)  The  SAR  image  was  recorded  at  17:44:29.  The  incidence 
angle  at  the  ship  track  was  47.6’.  The  appropriate  Bragg  wavelength  at  L-Band  (1.24  GHz)  is  16.5 
cm,  and  at  A-Band  (9.35  GHz)  is  2.2  cm. 

Attributes  leading  to  the  selection  of  JOWIP  1-1  were: 

a.  Relatively  low  uncertainties  in  crucial  parameters  such  as  the  pattern  propagation 
speed  and  ship  track  location; 

b.  Availability  of  a  nearly  contemporaneous  SAR  image; 

c.  Availability  of  good  slope  gauge  data  throughout  the  run; 

d.  Availability  of  video  and  real  aperture  radar  data  in  some  portion  of  the  run; 

e.  Simplicity  of  the  geometry:  the  near  alignment  (or  counteralignment)  of  wind, 
internal  wave  and  ship  direction; 

f.  A  distinctive  internal  wave  signal,  covering  a  range  of  magnitudes  (unlike  the  case 
of  JOWIP  12-3/4,  where  internal  waves  of  exceptionally  large  amplitude  were 
present). 

This  is  also  the  run  in  which  all  participating  agencies  have  expended  substantial  efforts  into 


mv 
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data  reduction. 

The  major  drawback  of  JOWIP  1-1  is  that  only  a  short  portion  of  real  aperture  radar  and 
video  camera  data  exist,  and  they  are  non  overlapping. 

i)  Calculated  Bragg  wave  modulation. 

Both  the  16  and  64  mode  smoothed  current  profiles  shown  in  Figure  2a  were  used  to  calculate 
the  spectral  modulation  for  waves  which  are  Bragg  resonant  with  the  ERIM  SAR  at  L-  and  X- bands  at 
the  look  angle  47.6*.  It  was  noticed  that  magnitude  and  phase  differed  significantly  in  the  calculated 
results  -  the  16  mode  profile  generally  yielding  a  smaller  perturbation.  In  this  paper,  all  calculated 
results  shown  were  obtained  using  the  64  mode  current  profile. 

The  calculated  results  at  kg  =  38.36  m-1  (X-band)  and  kg  =  289.22  m-1  (X-band)  are 
compared  in  Figure  3.  It  is  clear  that  the  predicted  X-band  modulation  is  significantly  smaller  than 
the  predicted  X-band  modulation:  the  range  of  variation  in  f  (i.e.,  fmaz  -  /  .  )  for  X-band  is  3.4  times 
larger  than  for  X-band. 

To  better  quantify  the  comparison,  we  define  the  averaged  squared  modulation  density 
D{X-,X0)  as 


«**•>  *  prin  J*  '"<•>  *• 


D{X\Xg)  is  a  measure  of  the  available  signal  averaged  over  the  interval  extending  from  x  —  X0  to 
x  —  X  and  gives  a  measure  of  the  potential  detectability  of  the  signal  within  the  interval.  In  these 
terms,  the  potential  detectability  for  the  current  record  extending  from  r=20mtox=  500  m  at  X- 
band  is  some  23  times  larger  than  at  X-band  (cf.  Figure  4).  (Note  that  the  main  current  structure  is 
confined  between  x  =  20  and  x  =  400.) 

ii)  Comparison  with  ERIM  SAR  Cut. 

Comparison  between  these  calculations  and  the  ERIM  SAR  cuts  are  shown  in  Figures  5  and  6 
for  X-  and  X-bands,  respectively.  What  is  compared  here  is  the  Bragg  model  prediction  F(k,x)/ Fef(k) 
with  the  observed  ratio  &(k,x)/&it(k),  where  i7t(  is  the  background  signal  as  approximated  b''  the 
mean  level  in  the  interval 

The  X-band  comparison  shows  that  the  calculated  modulations,  as  given  by  the  range  of 
variation  in  /,  are  larger  by  27%  than  the  modulations  in  the  ERIM  SAR  cut.  The  calculated  X-band 
results,  however,  underpredict  the  SAR  cut  by  a  factor  of  almost  3. 
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The  detectability  measure  D  defined  in  Equation  (13)  can  be  applied  to  the  case  of  cross 
section  data  by  replacing  /  in  its  integrand  with  the  relative  change  The  comparison  in  terms 

of  D  for  the  overlapping  interval  from  x  =  120  m  to  z  =  500  m  is  displayed  in  Figure  7.  The  overall 
SAR  realized  detectability  is  seen  to  be  comparable  for  Z-  and  A-band,  while  the  detectability 
predicted  by  the  Bragg  model  is  over  20  times  larger  for  Z-  than  for  A-band. 

From  these  comparisons,  one  is  drawn  to  the  conclusion  that  the  model  appears  to  perform 
reasonably  well  in  predicting  the  amplitude  of  the  Z-band  SAR  signal,  and  to  underpredict  the 
corresponding  A- band  SAR  result  by  a  factor  of  order  10.  (See  Table  I) 

What  is  calculated  here  is  the  maximum  signal  available  in  the  context  of  Bragg  scattering  to 
an  idealized  detector  with  no  loss.  A  more  meaningful  comparison,  for  the  purpose  of  validating  the 
model,  is  to  compare  the  results  to  surface  slope  modulation  data  and  real  aperture  radar  data. 


iii)  Comparison  with  Band-Passed  Slope  Gauge  Data. 

The  modulation  of  slope  in  the  Bragg  resonant  component  of  the  surface  wave  was  obtained  by 
band-pass  filtering  the  slope  gauge  data  at  the  appropriate  time-dependent  frequency.  This  time 
dependence  arises  from  the  need  to  account  for  the  Doppler  shift  induced  by  the  current  as  seen  by  the 
ship-mounted  slope  gauge.  This  is  equivalent  to  a  wavenumber  filter  over  a  range  [kg  —  &k-, 
kB  +  Ak+]  centered  at  the  Bragg  wavenumber  kg,  where  Ak. j_  is  about  15%  of  kg.  In  the  calculations, 
we  ignore  this  finite  bandwidth  and  assume  a  perfect  filter,  comparing  the  calculated  f(kg)  with  the 
band-passed  slope  results  for  Z-  and  A- band  Bragg-wave  modulations  (see  Kwoh  et  al,  this  issue). 

In  preparation  for  comparison  with  calculated  results,  a  21-point  running  average  (6.3  sec)  was 
performed  on  the  raw  data,  as  shown  in  Figure  8.  This  process  reduces  the  total  dynamic  range  of  the 
data  by  about  5  dB  in  this  case.  Comparison  between  the  calculated  results  and  the  band-passed  slope 
data  is  summarized  in  Table  II  and  displayed  in  Figures  9  and  10  for  Z-band  and  A- band,  respectively. 

It  can  be  seen  that  for  Z-band,  the  magnitude  of  the  calculated  result  is  in  good  agreement 
with  the  filtered  slope  data  record.  The  phase  of  the  modulations  is  not  in  total  agreement.  This  may 
be  due  to  the  inherent  uncertainties  associated  with  phase  speed  estimates,  unsteadiness  of  the  internal 
wave  features,  and  relative  location  of  the  current  and  slope  gauge  sensors.  The  agreement  should  be 
considered  acceptable. 


The  A-band  comparison  shows  that  the  calculations  underpredict  the  magnitude  of  the 
modulation.  The  underprediction  (4  db  predicted  dynamic  range  vs  25  db  observed),  is  much  more 
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serious  than  that  observed  when  the  calculations  were  compared  with  the  ERIM  SAR  cut  (4  dB 
predicted  vs  8  dB  observed).  In  terms  of  the  detectability  measure,  the  calculated  and  experimental 
D(500;20)  are  very  close  for  £-band.  For  Y-band,  D(500;20)  based  on  the  slope  gauge  measurements  is 
almost  24  times  larger  than  the  value  based  on  the  calculated  /. 

The  different  degree  of  agreement  between  the  calculations  and  the  data  (actually,  the 
difference  between  the  slope  gauge  and  SAR  measurements,  cf.  Figure  11)  is  consistent  with  the 
interpretation  that  the  band-passed  slope  data  represents  the  total  available  hydrodynamic  signal  under 
the  Bragg  assumption,  and  can  be  achieved  only  with  an  ideal  detector/processor  system.  It  also 
indicates  that  the  model  is  not  adequate  for  Y-band. 

iv)  Comparison  with  APL  Video  Data. 

The  APL  video  data  exists  for  a  short  range  of  JOWIP  1-1  and  is  compared  with  the 
calculated  results  in  Figure  12  for  the  L-band  case.  It  can  be  seen  that  the  comparison  is  good,  once 
again  confirming  that  at  least  for  JOWIP  1-1,  the  Z-band  results  are  well  predicted  by  the  calculations. 

v)  Comparison  with  TRW  Y-Band  Real  Aperture  Radar  Data. 

The  TRW  Y-band  real  aperture  radar  operates  at  9.13  GHz.  During  JOWIP  1-1,  it  was 
mounted  in  such  a  way  that  the  incidence  angle  was  40*.  Data  does  not  exist  over  the  entire  range  of 
JOWIP  1-1,  but  covers  a  time  period  starting  slightly  before  the  initiation  of  the  run  and  ending  before 
the  large  event  in  the  record. 

Comparison  with  the  ERIM  SAR  cut  (Figure  13)  indicates  that  the  Y-band  real  aperture  radar 
has  a  much  stronger  dynamic  variation  than  the  ERIM  Y-band  SAR  cut  (in  the  overlapping  region 
(from  x  =  120  m  to  x  =  280  m),  again  supporting  the  contention  that  there  may  be  more  signal 
available  than  that  detected  by  the  SAR. 

The  model  is  found  to  underpredict  the  RAR  measurements  by  a  factor  of  almost  10  in  the 
maximum  modulation  strength  found  in  the  range  x  =  20  m  to  z  =  280  m:  total  predicted  modulation 
of  35%  to  total  measured  modulation  of  about  270%  (cf.  Figure  14  and  Table  III). 

The  data  displayed  in  Figure  15  clearly  establish  a  one-to-one  relationship  between  real 
aperture  radar  backscatter  and  band-passed  slope  power.  In  particular,  the  values  for  D(280;20)  differ 
by  less  than  15%.  The  excellent  agreement  (observed  also  in  all  of  the  other  runs)  validates  the  band¬ 
pass  process.  Thus,  even  though  the  real  aperture  radar  data  is  lacking  over  a  significant  range  making 
direct  quantitative  comparison  with  calculations  difficult,  conclusions  on  real  aperture  radar 
backscatter  modulation  are  well  supported  by  comparison  using  the  band-passed  slope  data  in  this  case. 


vi)  Summary  for  JOWIP  1-1. 


Calculated  modulations  of  the  Bragg  resonant  component  using  TRW-HIC  with  the  Snyder  et 
al  (1981)  relaxation  rate  have  been  compared  with  experimental  data.  For  i-band  radar,  the  model 
appears  successful  in  predicting  the  measured  modulations,  with  perhaps  a  slight  underprediction.  For 
X-band  radar,  the  model  seriously  underpredicts  the  measurements. 

The  X-band  underprediction  appears  less  severe  in  comparisons  with  ERIM  SAR  cuts  than  in 
comparisons  with  band-passed  slope  data  and  real  aperture  microwave  data.  Part  of  this  may  be 
caused  by  the  fact  that  the  band-passed  slope  data  represents  the  maximum  available  modulation  due 
to  Bragg  backscatter,  and  should  only  be  achieved  by  an  ideal  detector/processor.  Since  the  model 
calculates  the  maximum  available  signal,  the  calculated  results  are  more  appropriately  compared  with 
the  band-passed  slope  data,  as  the  loss  in  converting  hydrodynamic  signature  to  SAR  processed  output 
is  not  included  in  the  model  used. 


b.  TRW  RUN  81B. 

This  run  corresponds  to  data  taken  in  the  time  interval  from  16:28:40  to  16:38:58  PDT  on 
August  4,  1983.  The  recorded  wind  speed  at  10  m  height  was  2.74  m/s  blowing  towards  123*T.  The 
ship  was  heading  3 ITT  at  an  approximately  constant  speed  of  1.2  m/s.  The  internal  waves  which 
induced  the  current  pattern  were  not  steady,  their  propagation  speed  ranging  from  about  0.2  m/s  to 
about  0.5  m/s,  with  an  average  speed  of  0.27  m/s  in  an  average  direction  128*T. 

A  reason  for  choosing  this  run  for  detailed  comparison  is  that,  although  no  concurrent  SAR 
data  is  available,  the  TRW  X-band  real  aperture  radar  data  is  complete. 

The  current  profile  used  in  the  calculation  and  the  wave  trajectories  of  interest  are  shown  in 
Figure  16. 

i)  Comparison  with  Band- Passed  Slope  Gauge  Data. 

The  band-passed  slope  data  for  L-  and  X-bands  is  displayed  in  Figure  17.  It  can  be  seen  that 
the  modulation  magnitude  is  comparable  for  the  two  cases.  In  fact,  the  X-band  modulation  seems  to 
be  slightly  stronger. 

The  comparison  between  the  L-band  calculated  result  and  data  is  shown  in  Figure  18.  The 
agreement  is  very  good,  even  though  there  is  a  moderate  underprediction  of  the  magnitude  of  the 
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modulation.  The  corresponding  comparison  for  the  A-band  case  (Figure  19)  shows  that  the  calculation 
grossly  underpredicts  the  modulation  by  more  than  one  order  of  magnitude.  Table  IV  summarizes  the 
relevant  information. 

ii)  Comparison  with  A'- Band  Real  Aperture  Radar  Data. 

Figure  20  compares  the  band-passed  slope  data  with  the  EE  real  aperture  radar  return.  In  this 
case,  the  dynamic  range  of  the  modulation  in  the  radar  return  is  somewhat  greater  than  that  in  the 
slope  data.  Calculated  modulation  and  radar  data  are  compared  in  Figure  21  and  Table  IV. 
Calculation  again  underpredicts  data  by  more  than  one  order  of  magnitude. 


iii)  Summary  of  TRW  RUN  81B. 

The  analysis  for  this  run  supports  the  general  conclusions  reached  for  JOWIP  1-1:  the  model 
appears  acceptable  at  I- band,  but  it  grossly  underpredicts  the  hydrodynamic  and  radar  cross  section 
modulations  at  A-band. 


4.  EFFECT  OF  LONG  WAVE  FIELD. 


The  wavelength  of  the  Bragg  resonant  surface  wave  component  is  about  2  cm  for  A- band  and 
17  cm  for  I-band.  One  of  the  reasons  proposed  for  explaining  the  poor  performance  of  the  model  at  A- 
band  is  that  these  short  Bragg  resonant  components  are  strongly  affected  by  longer  surface  waves. 
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A.  Formulation. 

Alpers,  Roes  and  Rufenach  (1981)  considered  three  mechanisms  by  which  a  deterministic  long 
wave  field  can  affect  the  statistical  properties  of  short  Bragg  waves  as  they  are  viewed  by  SAR.  These 
are,  the  tilt  mechanism,  in  which  the  long  waves  alter  the  local  angle  of  incidence  of  the  radar;  the 
hydrodynamic  mechanism,  in  which  short  wave  density  is  modulated  by  the  straining  action  of  the  long 
waves  (based  on  the  radiation  stress  mechanism);  and  the  velocity  bunching  mechanism,  in  which  the 
synthesized  radar  intensity  is  modulated  due  to  phase  modulation  in  the  Doppler  return  induced  by 
long  wave  orbital  velocity.  The  hydrodynamic  mechanism  affects  the  spectral  content  of  the  surface 
waves;  the  tilt  mechanism  affects  radar  backscatter  return;  and  the  velocity  bunching  mechanism 
affects  only  SAR  return. 

Alpers  et  al  showed  that  in  the  presence  of  a  deterministic  long  wave  field  described  by  the 
Fourier  spectrum  Z(k),  where  k  is  the  long  wave  wave-vector,  the  backscatter  due  to  Bragg  scattering 
can  be  expressed  as: 


Jk  dk  , 

Jo  J 


<r  =  «r>  1  + 


where  «r>  is  the  radar  backscatter  cross  section  when  the  long  wave  field  is  absent  or  averaged  out. 
The  integrand  is  defined  by 

(15)  Jk  -  R(i)  Z(k)  e +  c.c. 

■»  ^  ^ 
where  0(k)  is  the  long  wave  frequency,  Z(k)  is  related  to  the  surface  elevation  rjL(xrt)  of  the  long 

waves  through 


nL(*,t)  =  |  Z{k)  eik  UQt  dk  +  c.c.  , 


and  R(k)  is  a  complex  modulational  transfer  function  (MTF)  associated  with  the  long  wave  field  effect. 
The  MTF  R(k)  is  composed  of  three  terms  according  to  the  three  mechanisms  described  above: 
(17)  R  =  RTILT  +  RHYDRO  +  rvb . 

Detailed  expressions  for  R  can  be  found  in  Alpers  et  al  (1981)  or  Vesecky  and  Stewart  (1982). 

In  the  present  application,  all  of  the  waves  which  may  have  an  effect  on  the  short  Bragg  waves 
are  considered.  We  identify  «r>  with  the  Bragg  theory  cross  section  <r B  in  the  absence  of  long  waves, 
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and  combine  Equations  (14)  through  (17)  to  obtain 


(18)  <r  =  <rfl(l  +  A +  /,  +  /,) 


=  <rB  1  +  [  U  Jk(k)dk  +  f  R  Jk(k)dk  +  [  C  Jk(k)dk  . 

[_  Jo  >KV  iKR 


In  Equation  (18),  Ix  represents  the  contribution  from  waves  long  compared  with  the  current 
pattern  of  interest,  /2  represents  the  contribution  from  waves  longer  than  the  resolution  cell  width  but 
shorter  than  the  scale  of  the  current  pattern,  and  /3  represents  the  contribution  from  waves  longer  than 
the  “dividing  wavelength”  (empirically  defined  as  7  times  the  Bragg  wavelength,  where  7  is  order  1, 
usually  taken  to  be  between  2  and  6),  but  shorter  than  the  resolution  cell  width.  A  schematic  sketch  is 
shown  in  Figure  22. 


The  long  wave  field  is  taken  to  be  random  and  described  statistically.  The  RMS  of  a  taken 
over  spatial  scales  large  compared  with  the  resolution  cell  width  can  be  written  as: 


a  a  =  a , 


2  +  T*  +  J35 


Terms  linear  in  Ix  survive  the  averaging  process,  but  terms  linear  in  /2  or  I3  average  out  to 


""■W  JL 

zero.  The  long  wave  field  effects  are  represented  by  7j  and  /3  . 


In  order  to  compute  the  right  hand  side  of  Equation  (10),  we  assume  that  is  a  stationary 
random  process.  The  one-time,  two-point  correlation  function  p(£,f)  is  defined  as 


where  the  overbar  now  denotes  ensemble  average.  The  long  wave  power  spectrum  ♦^(i)  is  just  the 
«• 

Fourier  transform  of  p((): 


»*(*)  =  }  P(()  **  dl 


and  it  follows  that 


z(k)  r(k,)  =  *L(k)  sik-kj. 
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The  ensemble  averaging  can  now  be  identified  with  spatial  averaging  over  scales  long  compared 
with  the  resolution  cell  width  (and  the  Bragg  wavelength),  but  short  compared  with  the  spatial  scales 
of  the  surface  current  pattern  associated  with  subsurface  disturbances.  As  mentioned  earlier,  this 
averaging  process  leaves  unaltered,  but  averages  out  the  linear  contributions  of  /2  and  /3. 

From  the  definition  of  /2  and  /3,  we  obtain 


(Kr  (Kr  _ 

?  =  I  "  I  "i  Mh 

JKV  JKfj 

rKR 

=  |Jf(i)|3  #£(i,r)  dk 

1 K  rt 


where  the  dependence  on  *  is  made  explicit  to  emphasize  the  fact  that  the  long  wave  field  is  varying  on 
a  scale  large  compared  with  the  averaging  dimensions.  In  this  way,  the  modulation  due  to  the  current 

—  ■  i 

pattern  is  preserved.  Similarly,  /32  is  obtained  as: 

ft 

(24)  1?=  \  °  \RCk)\3*L(U)dl 

JKR 

The  modulated  RMS  radar  cross  section  that  corresponds  to  Bragg  backscattering  with  long 
wave  field  effects  included  is  therefore 

r  kc  -i 1/2 

(25)  a RMS  =  &  =  'A#  (!+/,)*  +  [  C  l*(*)|2  *j(i,i)  dK  , 

L  JKu 


with  Ii  given  by 


rKu 

h  =  f  Mh  dk. 

j  0 


In  Equation  (25)  the  term  involving  the  integral  represents  contributions  from  long  waves 
which  are  longer  than  the  dividing  wavelength,  but  shorter  than  the  current  pattern.  The  terms 
involving  7lt  on  the  other  hand,  represents  contributions  from  waves  comparable  to  or  longer  than  the 
current  pattern.  This  is  a  clutter  term,  since  it  produces  modulations  with  spatial  scales  similar  to  that 
of  the  signal  obtained  from  the  varying  current. 
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B.  Application  to  Present  Calculations. 

When  the  long  wave  field  effect  is  included,  the  baclcscatter  cross  section  in  the  absence  of  a 
current  becomes  (in  the  unidirectional  formulation) 

(27)  <rti  =  <r0(kB)  [1  + 

where  <r0(kB)  has  been  defined  in  Equation  (2)  and  tg(kB)  is  the  long  wave  field  effect  on  the  Bragg 
wave  component  in  the  equilibrium  state  (i.e.,  independent  of  z).  In  the  presence  of  a  current,  the  long 
wave  field  effect  L  is  dependent  on  z  and  is  given  by 

(28)  l(k,r)  =  f  C  | *(*,tf)|J  *l(K,z)  dK, 

*  o 

where  V L(K,x)  is  the  long  wave  spectrum  perturbed  by  the  current.  The  maximum  potential  relative 
change  in  radar  cross  section  becomes 


Comparison  between  predictions  by  the  Bragg  model  and  predictions  by  the  model  including 
long  wave  effects  involves  the  comparison  of  f{kB,z)  with  the  right  hand  side  of  Equation  (29),  which 
we  denote  as  the  dimensionless  modulation  m(k,z). 

In  the  cases  exhibited,  the  look  direction  and  wind  direction  are  generally  aligned,  so  that  the 
velocity  bunching  mechanism  am  be  neglected.  When  comparing  with  band-passed  slope  gauge  data, 
only  the  hydrodynamic  term  (RwyDfl0)  in  Equation  (17)  is  operative.  When  comparing  with  real 
aperture  radar  data,  both  the  hydrodynamic  and  tilt  (RTILT)  terms  are  included. 


C.  Results  of  Calculation  with  Long  Wave  Field  Effects. 

The  comparison  between  relative  cross  section  variations  predicted  by  Bragg  theory  with  and 
without  long  wave  effects  for  JOWIP  1-1  is  shown  in  Figures  23  and  24  for  L-  and  X- bands 
respectively.  The  long  wave  field  effect  increases  the  potential  detectability  Z?(500;20)  by  less  than  10% 
at  L-band,  whereas  the  increase  is  almost  300%  at  A- band  (cf.  Table  V). 

This  effect  is  also  observed  in  the  case  of  TRW  RUN  8 IB.  The  comparison  of  results  with  and 


without  long  wave  field  effects  for  Z-  and  X-bands  is  shown  in  Figures  25  and  26  respectively.  In  this 
case,  Xg  =  —729  m  and  X  =  —9  m.  The  enhancement  at  Z-band,  as  given  by  the  detectability 
measure,  is  less  than  5%,  but  the  enhancement  at  X-band  is  almost  200%.  Long  wave  field  effects  on 
the  predicted  modulations  are  summarized  in  Table  VI. 


Part  of  the  reason  that  the  long  wave  field  effect  is  stronger  at  X-band  than  at  Z-band  is  that 
more  waves  qualify  as  “long  waves”  affecting  the  short  X-band  Bragg  waves.  Another  reason  is  that 
the  “long  waves”  included  in  the  X-band  case  are  strongly  modulated. 


It  should  be  pointed  out  that  since  the  long  wave  field  effect  depends  on  the  shape  of  the 
surface  wave  spectrum,  it  is  also  a  function  of  wind  speed.  This  dependence  on  wind  speed  is 
complicated.  On  the  one  hand,  an  increase  in  wind  speed  shifts  the  spectrum  to  longer  waves  and 
make  more  “long  waves”  available.  On  the  other  hand,  an  increase  in  wind  speed  weakens  the 
modulation  of  all  waves.  Thus,  for  a  given  current  pattern,  there  should  be  a  wind  speed  for  which  the 
signal  is  maximum. 


Figure  27a  compares  the  calculated  X-band  results  (with  long  wave  field  effects)  to  X-band 
slope  data  for  JOWIP  1-1.  The  discrepancy  is  reduced  but  there  is  still  an  underprediction:  the 
Z)(500;20)  value  based  on  the  calculation  underpredicts  the  actual  detectability  found  in  the  slope 
gauge  signal  by  a  factor  of  6  (cf.  Figure  27b). 


The  comparison  with  the  ERIM  SAR  data  is  displayed  in  Figure  28,  which  shows  a  theoretical 
underprediction,  based  on  Z)(500;120),  by  a  factor  of  1.6. 


The  corresponding  comparisons  for  TRW  81B  are  shown  in  Figure  29  (against  band-passed 
slope  data)  and  Figure  30  (against  real  aperture  radar).  Based  on  D{— 729,— 9),  the  calculation 
including  hydrodynamic  long  wave  effects  underpredicts  the  slope  gauge  data  by  about  a  factor  of  50, 
down  from  around  93  predicted  by  the  Bragg  model.  The  calculation  including  both  hydrodynamic 
and  tilt  mechanism  underpreticts  the  detectability  of  the  RAR  data  in  the  same  range  by  a  factor  of  53 
( VV  data)  and  a  factor  of  49  {HH  data).  The  underpredictions  without  long  wave  effects  were  by 
factors  of  107  and  138,  respectively. 


In  summary,  long  wave  field  effects  do  not  alter  the  relatively  good  agreement  at  Z-band,  but 
significantly  reduce  the  discrepancy  at  X-band.  However,  even  with  the  long  wave  field  effects 
included,  the  underprediction  at  X-band  is  not  eliminated. 
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5.  SUMMARY  AND  CONCLUSIONS 

Model  calculations  of  surface  wave  modulation  by  surface  currents  have  been  compared  with 
field  measurements  made  with  SAR  (A-  and  i- bands),  real  aperture  radar  (A-band),  slope  gauge 
(band-passed  at  A-  and  L-band  Bragg  wave  frequencies),  and  to  a  limited  extent,  CCD  video  camera 
(processed  for  L-band  Bragg  wave  wavenumber).  The  following  observations  have  been  made: 

1.  Among  several  wind-relaxation  models  used  in  the  calculation,  the  model  proposed  by 
Snyder  ct  al  (1981)  yields  the  best  agreement  with  data.  The  Inoue  (1966)  model,  used  extensively  in 
the  past,  yields  modulations  an  order  of  magnitude  smaller  and  should  be  rejected. 

2.  A  theoretical  model  baaed  on  simple  Bragg  scattering  and  wind-relaxation  predicts 
modulations  measurable  at  A'-band  to  be  an  order  of  magnitude  or  more  smaller  than  those  measurable 
at  L-band.  This  is  in  disagreement  with  the  observations  (including  SAR,  RAR  and  slope  gauge  data) 
which  show  that  the  modulation  magnitudes  at  A-  and  L-bands  are  similar. 

3.  The  model  calculations  for  L-band  modulations  show  reasonable  agreement  with 
measurements.  Specifically,  the  calculations  agree  well  with  slope  gauge  data  and  CCD  video  camera 
data,  and  slightly  overpredict  the  SAR  modulations. 

4.  The  model  calculations  for  A-band  modulations  underpredict  the  measurements.  Based  on  a 
measure  of  overall  detectability,  the  underprediction  is  by  about  a  factor  of  10  when  compared  with 
SAR,  and  by  a  factor  of  20  to  40  when  compared  with  real  aperture  radar  and  slope  gauge  data 
smoothed  to  approximate  a  spatial  resolution  comparable  to  that  of  the  SAR. 

5.  A  new  model  which  incorporates  the  effects  of  modulated  long  surface  waves  on  Bragg  waves 
shows  that  the  long  wave  field  (LWF)  effect  enhances  the  theoretical  results  by  only  a  few  percent  at 
£-band,  but  by  100%  to  300%  at  X-band.  Therefore,  the  LWF  effect,  when  included,  reduces  the 
discrepancy  between  model  and  data  at  A-band  without  affecting  the  generally  good  agreement  at 
L-band. 

6.  Even  with  the  LWF  effects  included,  the  model  underpredicts  A-band  modulation  by  a 
factor  of  3  when  compared  with  SAR,  and  a  factor  of  10  to  30  when  compared  with  real  aperture  radar 
and  slope  gauge  data. 


It  should  be  emphasized  that  the  above  observations  are  based  on  a  limited  set  of  data  and 
comparisons.  Further  analysis  of  data  for  cases  with  weaker  surface  currents  and  a  wider  range  of 
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surface  conditions  is  underway  and  will  be  reported  subsequently. 
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TABLE  III:  JOWIP  1-1.  (z  range:  20  m  to  280  m) 

RAR  modulations  compared  to  Bragg  theory  predictions 


F/F" 

De  tec  tab 

dynamic 

at 

range 

r=280  m 

A-band  calculated 

0.854 

1.205 

1.5  db 

0.006 
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TABLE  V:  JOWIP  1-1.  (r  range:  20  m  to  500  m) 

Long  Wave  Field  effects  on  predicted  modulations. 
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dynamic 
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2.616 

26.9  db 

0.333 

+Hydro+Tilt(  VV) 

0.005 
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0.530 
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TABLE  VI:  TRW  81 B.  ( x  range:  -729  m  to  -9  m) 

Long  Wave  Field  effects  on  predicted  modulations. 
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A-BAND: 


Bragg 

0.511 

1.320 

4.1  db 

0.010 

+Hydro 

0.441 

1.509 

5.3  db 

0.019 

+Hydro-fTilt(  VV) 

0.435 

1.561 

5.6  db 

0.020 

+ Hydro+Tilt(  HB) 

0.402 

1.787 

6.5  db 

0.028 
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FIGURE  CAPTIONS. 

Figure  1.  a)  Relaxation  rate  /?  as  a  function  of  wavenumber  k  in  the  range  containing  Bragg 

resonant  waves  with  L-  and  A-band  radiation.  (Wind  speed  is  5  m/s  blowing  in 

-a 

the  k  direction). 

b)  A-band  results  for  the  conditions  of  JOWIP  1-1  using  various  relaxation  rate 
models. 

Figure  2.  JOWIP  1-1: 

a)  current  as  a  function  of  distance  along  the  ship’s  track.  Full  line:  64  mode  FFT 
smoothing.  Dotted  line:  16  mode  FFT  smoothing. 

b)  Selected  one-dimensional  ray  trajectories  in  k-z  space  for  the  current  (64  mode 
smoothing)  displayed  in  (a).  The  horizontal  lines  indicate  approximate  position  of 
Bragg  resonant  k  values  for  L~  and  A-band  radiation. 

Figure  3.  JOWIP  1-1:  Current  induced  modulation  of  the  relative  perturbation  spectral 

density  F(kg,z)/Fei(kg)  =  1  +  f(kg,x)  calculated  for  the  64  mode  current  of 
Figure  2. 

Full  line:  I- band  resonant  wave  {kg=  38.36  m~^). 

Dotted  line:  A-band  resonant  wave  (k^=289.22  m**). 

a)  linear  scale, 

b)  log  scale  in  db. 

Figure  4.  Potential  detectability  [cf.  Equation  (13)]  for  the  results  displayed  in  Figure  3. 

Figure  5.  JOWIP  1-1:  £-band.  Comparison  between  calculated  Bragg  wave  modulation 

F(kg,z)/ FCf(kg)  (dotted  line)  and  relative  cross  section  <r(r)/<re,  (full  line)  in  the 
ERIM  SAR  cut. 

a)  linear  plot, 

b)  log  plot  in  db. 

Figure  6.  JOWIP  1-1:  A-band.  Comparison  between  calculated  Bragg  wave  modulation 


F(kg,x)/Fei(kg)  (dotted  line)  and  relative  cross  section  <t(x)I<tb  (full  line)  in  the 
ERIM  SAR  cut. 

a)  linear  plot, 

b)  log  plot  in  db. 

Figure  7.  JOWIP  1-1:  SAR  achieved  detectability  (full  line)  compared  to  Bragg  model 

potential  detectability  (dotted  line). 

a)  £-band, 

b)  A-band. 


Figure  8.  JOWIP  1-1:  Band-passed  slope  gauge  data.  Dotted  line:  raw  data.  Full  line:  21 

point  average  (6.3  sec). 

a)  filtered  around  f-band, 

b)  filtered  around  A-band. 

Figure  9.  JOWIP  1-1:  L-Band.  Comparison  between  slope  gauge  and  calculated  Bragg  wave 

modulation.  Full  line:  relative  power  around  kg  =  38.36  m‘*  in  units  of  the 
average  power  for  the  interval.  Dotted  line:  Calculated  F(kB,x)/Feq(kB) 

a)  linear  plot, 

b)  log  plot  in  db. 

Figure  10.  JOWIP  1-1:  A-Band.  Same  as  Figure  9  caption,  but  iB=289.22  m'*. 

Figure  11.  JOWIP  1-1:  Comparison  between  normalized  slope  gauge  power  (full  line)  and 

normalized  SAR  cross  section  (dotted  line). 

a)  f-band, 

b)  A-band. 


Figure  12.  JOWIP  1-1:  Normalized  spectrum  F(kg,x)/ Fei(kg)  at  £-band. 

Full  line:  APL  video  data  results, 

Dotted  line:  predicted  by  Bragg  model. 


Figure  13.  JOWIP  1-1:  A-band.  Comparison  between  TRW’s  RAR  (full  line)  and  ERIM’s 

SAR  (dotted  line)  cross  sections. 


Figure  14. 


JOWIP  1-1:  A-band.  Comparison  between  normalized  RAR  cross  section  (full 
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line)  and  Bragg  model  normalized  spectral  density  (dotted  line). 

a)  linear  plot, 

b)  log  plot  in  db. 


Figure  15. 


Figure  16. 


Figure  17. 


Figure  18. 


Figure  18. 


Figure  20. 


Figure  21. 


Figure  22. 


JOWIP  1-1:  X-band.  Comparison  between  normalized  RAR  cross  section  (full 
line)  and  normalized  slope  gauge  power  (dotted  line). 

a)  linear  plot. 

b)  log  plot  in  db. 

TRW  81B: 

a)  Surface  current  profile. 

b)  Selected  ray  trajectories  in  k-x  space  for  the  current  shown  in  (a). 

TRW  81B:  Band  passed  slope  gauge  data.  Full  line:  £-band;  dotted  line:  X-band. 

a)  linear  plot, 

b)  log  plot  in  db. 

TRW  81B:  L-band.  Comparison  between  band-passed  slope  gauge  data  (full  line) 
and  Bragg  Model  Calculation  (dotted  line). 

a)  linear  plot, 

b)  log  plot  in  db. 

TRW  81B:  X-band.  Comparison  between  band-passed  slope  gauge  data  (full  line) 
and  Bragg  Model  Calculation  (dotted  line). 

a)  linear  plot, 

b)  log  plot  in  db. 

TRW  81B:  X-band.  Comparison  between  band-passed  slope  gauge  data  (full  line) 
and  TRW’s  RAR  normalized  HE  cross  section. 

TRW  81B:  X-band.  Comparison  between  TRW’s  RAR  Rff  cross  section  (full  line) 
and  Bragg  Model  Calculation  (dotted  line). 

Schematic  representation  of  long  wave  contributions  to  the  spectrum  modulation 

at  kg. 
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Figure  23. 


Figure  24. 


Figure  25. 


Figure  26. 


Figure  27. 


Figure  28. 


Figure  29. 


] 

JOWIP  1-1:  L-band.  Comparison  between  computational  results  with  (full  line) 
and  without  (dotted  line)  long  wave  field  effects. 

a)  predicted  relative  cross  section, 

b)  predicted  potential  detectability. 

JOWIP  1-1:  Af-band.  Comparison  between  computational  results  with  (full  line) 
and  without  (dotted  line)  long  wave  field  effects. 

a)  predicted  relative  cross  section, 

b)  predicted  potential  detectability. 

TRW  81B:  L-band.  Comparison  between  computational  results  with  (full  line) 
and  without  (dotted  line)  long  wave  field  effects. 

a)  predicted  relative  variation  in  cross  section, 

b)  predicted  potential  detectability. 

TRW  81B:  A"- band.  Comparison  between  computational  results  with  (full  line) 
and  without  (dotted  line)  long  wave  field  effects. 

a)  predicted  relative  variation  in  cross  section, 

b)  predicted  potential  detectability. 

JOWIP  1-1:  Af-band.  Comparison  between  band-passed  slope  gauge  data  (full 
line)  and  LWF  model  calculations  (dotted  line). 

a)  normalized  power, 

b)  detectability. 

JOWIP  1-1:  Af-band.  Comparison  between  ERIM’s  SAR  data  (full  line)  and  LWF 
model  calculations  (dotted  line). 

a)  relative  cross  section, 

b)  detectability. 

TRW  81B:  Af-band.  Comparison  between  ERIM’s  SAR  data  (full  line)  and  LWF 
model  calculations  (dotted  line). 

a)  relative  cross  section, 

b)  detectability. 


1 


Figure  30. 


TRW  81B:  Af-band.  Comparison  between  TRW’s  BE  RAR  data  (full  line)  and 
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LWF  model  calculations  (dotted  line). 

a)  relative  cross  section 

b)  detectability. 


Video  (L  Band) 
Bragg  Model 


FIGURE  24 


Band-passed  Slope  Gauge  Data 
LWF  Model  (hydro  only) 


